Backlund transformations for the sl(2) 

Gaudin magnet 



A N W Hone,t[] V B Kuznetsov§g and O Ragniscofg 

^ ' fDepartment of Mathematics, University of Kent, Canterbury CT2 7NF, UK 

(3 . §Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT, UK 

^ ! JDipartimento di Fisica 'E. Amaldi', Universita degh Studi 'Roma Tre', Roma, Itaha 

O 
D 

Q 

in 

I ■ 

Abstract 

! Elementary, one- and two-point, Backlund transformations are constructed for the generic 
case of the sl(2) Gaudin magnet. The spectrality property is used to construct these explicitly 
given, Poisson integrable maps which are time-discretizations of the continuous flows with 
^ ■ any Hamiltonian from the spectral curve of the 2x2 Lax matrix. 
> ' 

O 
l> 

o 
o 
o 

> 
X 



^E-mail: A.N.W.Hone@ukc.ac.uk 
^E-mail: vadimOainsta. leeds.ac.uk 
•^E-mail: ragniscoOf is .uniromaS . it 

1 



1. Introduction 



Backlund transformations (BT's) are an essential tool used to generate new solutions out of 
given solutions to integrable equations. This is by now a well-developed area, with elegant 
BT's having been found and studied for almost all integrable hierarchies, see 0. 

The theory of BT's for evolution equations had entered the subject of finite-dimensional 
integrability through the discretization of time variable(s). One of the most important and 
earliest accounts on this subject is the papers by Veselov where integrable Lagrange 
correspondences were introduced as discrete-time analogs of integrable continuous flows. 
Veselov clarified geometric meaning of these correspondences as finite shifts on Jacobians 
and gave several important examples. Reader is referred to an extensive literature which 
has appeared since then, see for instance [§, ||, H, and references therein. 

In the present paper, following the approach of we look at BT's for finite-dimensional 
(Liouville) integrable systems as special canonical transformations, thereby taking a Hamil- 
tonian point of view. We introduce and study several new properties of BT's which appear 
to be very natural in such approach. 

Backlund transformations for finite-dimensional integrable systems are defined in this 
paper as symplectic, or more generally Poisson, integrable maps which are explicit maps 
(rather than implicit multi- valued correspondences of 0) and which can be viewed as time 
discretizations of particular continuous flows. The most characteristic features of such maps 
are: i) BT's preserve the same set of integrals of motion as does the continuous flow which 
they discretize, ii) they depend on a (Backlund) parameter A that specifies the corresponding 
shift on a Jacobian or on a generalized Jacobian 0, iii) a spectrality property holds with 
respect to A and to the 'conjugate' variable /z, which means that the point (A, /i) belongs to 
the spectral curve p, |10 



Because of the above properties, the constructed BT's are suitable as explicit (symplec- 
tic) geometric integrators. Namely, explicitness makes these maps to be pure iterative, while 
the importance of the parameter A is that it allows an adjustable discrete time step. The 
spectrality property is strongly related to the simplecticness of the map. Finally, numer- 
ical integrators which exactly preserve the level set of integrals and at the same time are 
symplectic proved to be impossible to find for generic Hamiltonian dynamics but for 
integrable flows they do exist and so are in demand. 

In this paper we consider a generic (diagonal) case of the sl(2) XXX Gaudin magnet 
which is an algebraic completely integrable system. We study the problem of constructing 
elementary (one- and two-point) Backlund transformations for this system. In Section 3 
we construct an elementary (one-point) BT which gives an exact discretization of a specific 
continuous flow. By making a two-point composite map in Section 7 we are then able to 
discretize any of the independent commuting flows with the Hamiltonians from the spectral 
curve of the 2x2 Lax matrix. 



2. Gaudin magnet 

The sl(2) Gaudin magnet is derived from the Lax matrix 



A{u) = a + Y^ 



B{u) = J:^, CH = E^- (2.2) 

Pi u - aj p^ u - 

Local variables in this model are generators of the direct sum of n sl(2) spins, s^, s^, 
j' = 1, . . . , n, with the following Poisson brackets: 
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We denote Casimir operators (spins) as sf 



^0 



(2.3) 



(2.4) 



Fixing Casimirs we go to a symplectic leaf where the Poisson bracket is non-degenerate, 
so that the symplectic manifold is a collection of n spheres. 

Let us also introduce the total spin J which will be used later, as follows: 



The Lax matrix ( |2.1|) satisfies the linear r-matrix Poisson algebra, 
{Li(m), L2{v)} = [r{u - v), Li{u) + L2{v)], 
with the permutation matrix as the r-matrix 



r{u — v) 



u — V 



/I 0\ 

10 

10 
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Here we use standard notations Li and L2 for tensor products: 



Li{u) = Liu) 



1 
1 



1 
1 



Uv) 



(2.5) 



(2.6) 



(2.7) 



(2.^ 



Equation ( [^.6D is equivalent to the following Poisson brackets for the rational functions A{u) 
B{u), and C{uy. 



{A{u),A{v)} = {B{u) 


,B{v)} = {C{u),C{v)} = 0, 


(2.9) 


{A{u),B{v)} = 


-^{B{v)-B{u)), 

U — V 


(2.10) 


{A{u),C{v)} = 


-^{C{v)-C{u)), 

u — V 


(2.11) 


{C{u),B{v)} = 


^{A{v)-A{u)). 
u — V 


(2.12) 



The spectral curve F, 



F: det(L(M) -t;) = 0, 



:2.13) 



3 



is a hyperelliptic, genus n — 1 curve, 



= A\u) + B{u)C{u) =a^ + Y. 



with the Hamiltonians (integrals of motion) Hj of the form 



(2.14) 



aj — ttk 



(2.15) 



These are integrals of motion, or Hamiltonians, of the sl(2) Gaudin magnet, which are 
Poisson commuting: 

{H^,Hk} = 0, j,fc = l,...,n. (2.16) 
Notice that there is one linear integral: 



We can bring the curve F into the canonical form by scaling the variable v ^ v: 



[2.17) 



V = V Y\_[u — a 



'J J- 



(2.18) 



The equation of the curve becomes 



a 



+ E 



+ 



^ \u — aj {u — aj) 



X{{u-aj 



2n-l 



,2n-2 



2n- 



(2.19) 



When a = the genus of the curve drops to n — 2, because /i = in such case. The Gaudin 
magnet then becomes sl(2)-invariant: apart from integrals ( |2.15 ) all three components of the 
total spin J are integrals too. 



« = 0: {H„Jk} = 0, j = l,...,n, k = 1,2,3. 



(2.20) 



We will not consider this case, but concentrate on the generic case of a 7^ when there is only 
one linear integral /i = 2a(J3 — aJ2]=i o^j)- The latter case is called generic (diagonal) case 
of the sl(2) XXX Gaudin magnet. It is known that all its flows are linearized on generalized 
Jacobian of the hyperelliptic curve ( |2.19| ), see references in 0. 



3. One-point basic map 

The sl(2) Gaudin magnet with the 2x2 Lax matrix ( ^.1| ) is within the class of systems that 
was considered recently in 0, namely it belongs to the (even) case of the generalized Jaco- 
bian. Hence, its Backlund transformations can be extracted from that paper. However, we 
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want to present here an independent derivation of those BT's as well as to give more detailed 
exposition of their various properties. The reader is referred to also for explanation of 
geometric meaning of Backlund transformations. 

A Backlund transformation should act on the Lax matrix as a similarity transform: 

L{u) ^ M{u)L{u)M{u)-^ Vw, (3.1) 

with some non-degenerate 2x2 matrix M(m), simply because a BT should preserve the 
spectrum of L{u). 

Let us introduce new (~-) notations for the up-dated variables 

_ n~3 _ n ~ - _ n~ + 

I(n) = a + E^^, 5(n) = E^^, C{u) = Y.^^, (3.3) 

j=l u - j=l u - «j j=l u - «i 

{s/, Sk^} = T^SjkSk^, {sj^, Sk'} = -2\5jkSk. (3.4) 
We are looking for a Poisson map that intertwines two Lax matrices L{u) and L{u): 

M{u)L{u) = L{u)M{u) Vu. (3.5) 

Because spins Sj, j = 1, . . . ,n, appear as coefficients of the curve, they are not changed by 
the map, i.e. Sj = Sj. Hence, we can talk about a symplectomorphism [sj = const) instead 
of a Poisson map. 

Now we should choose an ansatz for the dependence of the matrix M{u) on the spectral 
parameter u. Let us fix the simplest case of a linear function: 

M{u) = Miu + Mq. (3.6) 

Taking limit m — >■ oo in (|3.5|) we conclude that Mi must be diagonal. Moreover, the most 
elementary (one-point) BT should correspond to the case when det M{u) has only one zero 
M = A, which will lead to having only one Backlund parameter (cf. the spectrality property 
in H, |l^). So, we should choose either 

M.-d o. M.^f« °V ,3.7) 



OJ ' vo 1 

We will consider the first the second one will produce a similar BT (cf. moving 

in discrete time in positive and in negative direction). Finally, we arrive at the following 
parameterisation of the unknown matrix M{u): 

MM=("-^ + «/^ P), (3.8) 

with the c-number determinant 

detM(M) = 7(m- A). (3.9) 
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Here the variables p and q are indeterminate dynamical variables, but A and 7 are c-number 
Backlund parameters^. 

Comparing asymptotics in u in both sides of ( p.5|) we readily get 

(3.10) 

If we want an explicit single- valued map from L{u) to L{u) 

L{u) ^ L{u) = M{u)L{u)M-\u), (3.11) 
then we must express M{u), and therefore p and q, in terms of the old variables, i.e. the 



entries of L{u) only. There is however a problem, since from (|3.1CI| ) we have only the expres- 
sion for the p, but the variable q is given in terms of new, and therefore unknown variable 
J+. To overcome this difficulty we will use an extra piece of data, namely the spectrality. 
Apart from equations ( |3.5| ) that our map satisfies, it will be parameterised by the point P 
on the curve F, 

P={X,IJ,)eT. (3.12) 

Notice that there are two points on the curve T, P = (A, +/i) and Q = (A, — /i), corresponding 
to a fixed A and sitting one above the other because of the hyperelliptic involution: 

(A,/i)Gr: det(L(A) - ;u) = ^ /i^ + det(L(A)) = 0. (3.13) 

Because det(M(A)) = the matrix M(A) has a one- dimensional kernel 

Mm={'f ^)f! = 0, fi=(:J. (3.14) 

The equality 

M{\)L{\)n = L{\)M{\)n = (3.15) 

implies that L{X)Q ~ Q, so that Q is an eigenvector of L{X). Let us fix the corresponding 
point of the spectrum as P = (A,/i): 



fAi\)-fx 5(A) \/7 
I C(A) -A{X)-fij{-q 



0. (3.16) 



This gives us the formula for the variable q: 



The two last expressions are equivalent since (A,/i) G F. 

Now, the formulas ( |3.11| ), (|3.8|), (|3.10| ), and ( p.l7| ), give a one-point Poisson integrable 
map (= one-point BT) from L{u) to L{u). The map is parameterised by one point (A,/i) on 
the spectral curve F (and by an extra parameter 7). 

Explicitly, it reads 

~ -f{u-X + 2pq/-f)A{u)-q{u-X+pq/-f)B{u)+p-fC{u) 

A[u) = , (3.18) 

7(m - A) 



^M(u) comes from a L-operator of the quadratic r-matrix algebra whose det is a Casimir, cf. E 



~ {u-X + pq/-tfB{u)-2p{u-\ + pq/-i)A{u)-p^C{u) 

= ' (3-19) 

7(m - A) 

Equating residues 0.1 u — aj in both sides of the above equations, we obtain the map in 
terms of the local spin variables: 

~ 3 _ 7(% - A + 2pg/7)gj _ g(a^. _ A + pq/^)sj + P7g+ 

~ - _ (% - A + pqhfs- - 2p{aj - A + pg/7)gf - P^s+ 
Sj — / w 1 \o.ll) 

li^j - A) 

7^4 + 297^,^ - q^s- 
- 7(a,- - A) ■ ^"^-^"^^ 

Recall that a, and Sj, j — 1, . . . ,n, are parameters of the model; 7 and A are parameters 
of the map; and variables p and q are as follows: 

J_ ^(A)-// C(X) , 

= ^ + ,^+,- + ^ ^^^^^^ ^^^^^^ 



4. BT as a discrete-time map 



In this Section we will show that the BT constructed above can be seen as a time-discreti- 
zation of a specific Hamiltonian flow where the parameter A plays a role of inversion of the 
time step. 

Consider the limit A ^ 00 then 



// = a + oQ). (4.1) 



Assume that 



7 = -A + 7o + oQ), (4.2) 
then we have the following expansion for the matrix M{u): 

M{u) = -A (1 - ^Mo(«)) + O (i) . (4.3) 
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The equation of the map, M{u)L{u) = L{u)M{u), turns in this hmit into the Lax pair of a 
continuous flow: 

L{u) = [L{u),Mo{u)], M,{u)={^j'J^ -t^To)' ^^'^^ 

where 1/(2A) is a time-step and L{u) = hm^^oo 2A(I/(m) — L{u)) is the time-derivative. 
The flow (\iA) is a Hamiltonian flow 



L{u) = {H,L{u)}, (4.5) 

with the Hamiltonian function H as 

H= - I J+J_ + 2af](a,-7o)4) . (4.6) 



Therefore, the constructed Backlund transformation is a two-parameter (A, 7) time-discreti- 
zation of this continuous flow. 



5. Symplecticity 

In this Section we give a simple proof of symplecticity of the constructed map by finding an 
explicit generating function of the corresponding canonical transformation from the old to 
new variables. 

First, because the spin variables (Casimirs) do not change. 



(4) +4^7 = (■'/) (5.1) 

we can exclude the variables and Sj^ , j = 1, . . . ,n, 

4 = ^^, sf = ^^^, (5.2) 

expressing everything in terms of 2n 'canonical' variables (s|,s~)"^i and (sj^, Sj"*")"^^ with 
the following Poisson brackets: 

} = i^jfcSfc , {S'i^ Sk'^} = -i6jkSk^. (5.3) 

We want to represent our Backlund transformation as a canonical transformation defined by 
the generating function F(§'"'"|s~) = Fa,^(s "'"|s~) such that 



3 ._dF{s+\s-) ^3 .^,dF{s+\s-j , ^, 



Notice that we have chosen the arguments of the generating function as {sj~^\Sj )]=i- Because 
the symplecticity property does not depend on the choice of the arguments of its generating 
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function, these arguments are fixed in order to get a simpler expression for the function F\^^ 
(recall that the variables p and q ( p.lO|) depend exactly on these variables). 
Rewrite now the equations of the map ( p.21|) -( ^.23| ) in the form 



74-^^7) + 7(% - A)s/s7 - = 0, (5.5) 



2a 



(73^/ - § + 7(% - A)?/s7 - ^^s] = 0. (5.6) 



Resolving them with respect to and Sj^, we obtain 



a, — A 

7 

It is now easy to check that the function 



^'^3 = ^? - ^J^^i ' J = 1, . . . , n. (5.8) 



satisfies equations (|5.4|). Thereby symplecticity of the map is proven. 



6. Spectrality 

The map depends on two parameters, A and 7. Let us first concentrate on its A-dependence. 

Spectrality, which was introduced in is an interesting property of Backlund transfor- 
mations. It usually holds for any Backlund transformation which has a parameter. Techni- 
cally, it means that the components of the point P = (A,/i) G F which parameterises the 
map are conjugated variables, namely: 

^^ aFU?1.-) ^ (6.1) 

To prove this formula, use ( p.lO| ) and ( |3.17| ) to find the formula for the 



/i = A(A)-^5(A). (6.2) 



Now, with the help of (|5.7|) and ( |5.9| ) we easily check the needed formula for the spectrality 
property ( |6.1|) . 



A new, comparing to [Q, observation is that there is also an analogous property with 
respect to the parameter 7, only now it is the integral J3 that plays the role of the conjugated 
variable: _ 

h = -7 • (6.3) 
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The proof is very simple, once we notice that ( |5.7] ) entails 



Concluding this Section we want to remark that because of the second 'spectrality' property 
( |6.3|) , which was somehow built into the Backlund transformation from the very beginning, 
one could recover the generating function of the map just taking one integral. 



f — — I d7 + const, (6.5) 
without needing to solve the system of 2n differential equations ( |5.7| ). 



7. Inverse map and a two-point map 

In this Section we will first construct the inverse Backlund transformation and then use it 
to derive the two-point Backlund transformation which will be a composition of the direct 
map parameterised by the point Pi = (Ai,/ii) G T and the inverse map parameterised by 
the point Q2 = (A2, -Ai2) G T. 

7.1 The inversion of the Backlund transformation 

Let us call the direct map by Bp. The inverse map acts from L{u) to L{u). We can rewrite 
the equations for the Bp in the inverse form 

M\u)L{u)=L{u)M'{u), M\u) = ( u-\lpqh) " ^^"^^ 

To define the inverse map we must find expressions for the co-factor matrix M^(-u), or for 
the variables p and g, in terms of ~- variables, i.e. in terms of the entries of L{u). We have 
already the expressions ( |3.1CI| ), 



which define q. To obtain the formula for the variable p we will use again the spectrality 
property. The matrix M^(A) has a one-dimensional kernel f2. 

The main difference comparing to the formulas of the direct map is that the inverse map 
will be parametrized by the point Q = (A, —fi) G F, not the P = (A,/i) G T. Therefore, fl is 
an eigenvector of the matrix L{u) with the eigenvalue Q = (A, — yu): 



{AiX)+fi 5(A) \fp 
V C(A) -A(A)+/iA7 



0. (7.4) 
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This gives us the needed formula for the variable p, 

p = 'y—^ = -7 ^ ^ ^ ■ (7.5) 

To prove that this does indeed give the inverse map, we have to show that the two formulas, 
( |3.17D and (|7.5|), define in fact the same variable /i: 

^ = A{\)-^C{\) = A{\)-^ B{\). (7.6) 

7 7 

It is easy to see that this equation is the (ll)-element of the matrix identity: 

M^(A)Z(A) = L(A)M^(A). (7.7) 

We will denote as Bq the map which is inverse to the map Bp. Generally speaking, we have 
constructed 4 different maps. Bp, Bq, Bq, and Bp, with two pairs of maps which are inverse 
to each other: 

BqoBp = BpoBq = BpoBq = BqoBp = Id. (7.8) 
7.2 The two-point map Bp^ q^ 

We now construct a composite map which is a product of the map Bp^ = B(^Xl,^J.l) and 

Bp,,q, = Bq,oBp,. (7.9) 

The second parameter of the basic map, namely the 7, is taken the same in both maps, so 
7i = 72- Obviously, when Ai = A2 (and /xi = /i2) this composite map will turn into an 
identity map. 

The first map Bp^ reads as follows: 

M,iu)Liu) = Liu)M,iu), Mi(t.) = ("-^^+^^^^/^ ) , (7.10) 

where the formulas for the variables pi and qi are 

^J_^ I(Ai)-/ii _ ^ ^(Ai) 

2a ^ C(Ai) ^ A(Ai)+/ii ' ^ ' ' 

J+ A(Ai)-^i CjX,) 
'' = 2-a='^{X^ = -'A(X^,- ^^-^2) 

The second map Bq^ reads as follows: 

M.{u)Uu) = Uu)M,{u), M,[u) = [2^^ u-X~+p,<i,h)^ ^'-''^ 
where the formulas for the variables p2 and q2 are 

^ ~J-^ AjXj) -112 _ ^ B{\2) 
2a ^ C(A2) ^/l(A2)+/X2' 
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B{\2) Ai\2)+fX2 

Notice that qi is equal to q2, hence we omit the sub-index, qi = q2 = q- 
The composite map Bp^ q^ acts from L{u) to L{u), 

M{u)L{u) = L{u)M{u), (7.16) 

M(u) = - AMu)AUu) = ^ , „^ _ _ „ _ _ , _ j . (7.17) 

In order to get rid of the intermediate ~-variables, we will use the spectrality property 
with respect to two points, Pi = (Ai,/ii) and Q2 = (A2, — /i2)- Obviously, both spectralities 
are still valid after composing the maps. For the point Pi we get the following equations: 

M(Ai)fii=0, ni=(^2q)^ L{Xi)ni = fii^i 
P(Ai) A{\i) + Hi 



M\Xi)Qi = 0, ^1 = ( A, - A^- Mpi - P2) ) ' ^ 



7(A,-Ai)g(Ai) li>^i->^2) (i(Ai)-/ii 

Pl-P2 = ^^ X ~ = \ ~ • 

7 (^A(Ai) + /ii j - gP(Ai) g (^A(Ai) - /^i j + 7C^(Ai) 
For the point Q2 we get the second set of equations: 

M{X2)Q2 = 0, Q,= (^^^_ _ ) ' ^(^2)^^2 = -/i2^^2 ^ 

„ „ _ 7(A2 - Ai)P(A2) _ 7(Ai - A2) (A(A2) - (^2) , . 

7(^(A2)+/i2)-gP(A2) g(A(A2)-/i2)+7C(A2)' ^' ^ 



M^(A2)^^2 = 0, n2=i^^^j, L{X2)n2 = fl2^2 

q^,m^^-,J^. (7.21) 

B{X2) A{X2)+fl2 

Equations ( [7.18| ) and ( |7.21| ) are already known to us (cf. ( |7.12| ) and (|7.15|) ). The formulas 
( f7l9D and d7!20| ) for the variable Pi —P2 are new. They are equivalent to the formulas ( [7. 11]) 

and ( [7. 141) expressed in terms of entries of L{u) and L{u). 

Concluding, we have constructed a two-point Backlund transformation which is factorised 
to two one-point Backlund transformations and which is explicitly given, together with its 
inverse, by the formulas: 

M(..)L(.) = £(..)M(..), M(u) = ( J-^''^l _ _ ) , (7.22) 

12 



where 



■ S(Ai)(A(A2)+/^2)-5(A2)(A(Ai)-/.i) ' ' 



(Ai-A2)i?(Ai)(A(A2)-/x2) 
(A(Ai) - /ii)(A(A2) - //2) + fi(Ai)C(A2) 

(A2-Ai)i?(A2)(A(Ai)+/xi) 
(A(Ai) + /xi)(A(A2) + 1^2) + 5(A2)C(Ai) 

(Ai-A2)(A(Ai)+^i)(A(A2)-/i2) 
(A(Ai) + /ii)C(A2) - (A(A2) - /x2)C(Ai) 

(A2-Ai)l(A2)^(Ai) 

^(A2) (l(Ai) + /ii) - ^(Ai) (1(A2) - /i2 

(Ai-A2)1(A2) (I(Ai)-/xi) 



(7.25) 
(7.26) 
(7.27) 

(7.28) 
(7.29) 
(7.30) 
(7.31) 

The above formulas give several equivalent expressions for the variables x and X since the 
points {Xi, III) and (A2, — A*2) belong to the spectral curve F, i.e. are bound by the following 
relations: 

= A^iXi) + i?(Ai)C(Ai), IJ.I = A\X2) + i?(A2)C(A2), (7.32) 
f,l = l'(Ai) + ^(Ai)§(Ai), f,l = l'(A2) + ^(A2)^(A2). (7.33) 

7.3 Two-point map as a discrete-time map 

We will see in this Section that the two-point map constructed above is a one-parameter, 
Ai, timc-discretization of a family of flows parameterised by the point Q2 — (A2, — /^2), with 
the difference Ai — A2 playing the role of the time-step. 
Indeed, consider the limit Ai — > A2, 

Ai = A2 + £, £ ^ 0. (7.34) 



^(A2) - 1^2) [A(Xi) - /^i j + S(A2)C(Ai) 

(A2-Ai)^(Ai) (1(A2)+/X2) 
MX2) + /i2) (i(Ai) + /xi) + ^(Ai)^(A2) 

(Ai - A2) (1(A2) + (l(Ai) - //i) 
1(A2) + /X2) ^(Ai) - (l(Ai) - /xi) ^(A2) 
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It is easy to see from the formulas of the previous subsection that 

A{X2)-fl2 C(A2 



X = xo + 0{e), 



Xq 



and 



X = eXo + 0{e^ 



B{\2) A{\2)+^i2 



2/i2 



Then we derive that the matrix M{u) has the following asymptotics: 

^' ^ 2/X2(m-A2) V C{\2) -A(A2)+/i2^ 

If we now define the time-derivative L{u) as 



0(e' 



(7.35) 



^7.36) 



(7.37) 



L(«)=lim^M^i^. 

£^0 F 



(7.38) 



then in the limit we obtain from the equation of the map, M{u)L{u) = L{u)M{u), the Lax 
equation for a corresponding continuous flow that our Backlund transformation discretizes, 
namely: 

LiX2) 



Liu) 



This is a Hamiltonian flow with /i2, 



Liu) 



2j22{u - A2) 



(7.39) 



/^2 = V^2(A2)+5(A2)C(A2) 

as the Hamiltonian function, 



\ 



a' 



+ E 



+ 



i=i \^2- a-j (A2 - a-j) 



L{u) = -i{/i2,L(M)}. 



(7.40) 



This means that the two-point map discretizes a one-parameter family of flows. Having 
chosen the parameter A2 to be equal to any of the poles of the Lax matrix (parameters of 
the model) aj, j = 1, . . . ,n, the map leads to n different maps, each discretizing the flow 



with the corresponding Hamiltonian Hj, j = ] 

A2 = CLj + £, 

Then we have 



, n. Indeed, take the limit Ao 



/i2 = ^ + ^ + 0(.), 
£ ZSj 

and in this limit the Lax equation ( |7.39| )-( |7.4CI| ) turns into 

1 



L{u) = ~—{H,,Liu)} 



Lu) 



2sj[u — ttj^ 



4 -4 



(7.41) 
(7.42) 

(7.43) 



Let us denote a collection of these maps by {Bp^}^^^. The map Bp^ discretizes the 
flow governed by the Hamiltonian with Ai — playing the role of the discrete time-step 
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parameter. The map (and its inverse) is defined by the two-point matrix M{u) ( [7. 221 ) with 
the following expressions for the variables x and X: 



. = (7.44) 

^ (Qfc - Ai)^(Ai)gfc (gfc - Ai)^(Ai)gfc 

5(Ai)(4 + ^.)-5.-(^(Ai)-/ii) ~s,- (A{X,) + f,,) - hx,) - s,] ' 



All these maps are explicit Poisson maps, preserving Hamiltonians and having the spec- 
trality property with respect to the pair of variables (Ai,/ii). 



8. Concluding remarks 

One of the very important branches of the theory of finite-dimensional integrable systems 
is the area of discrete-time integrable systems. The interest to this area was revived in 
the beginning of 90 's by Veselov in the series of works (see |^). He defined integrable 
Lagrange correspondences as discrete-time analogs of integrable continuous flows, clarified 
their geometric meaning as finite shifts on Jacobians and gave several important examples. 
Since then the subject has got a boost and has been developed further by many authors. It 
has not been our intention to give a review of many important recent contributions made 
to the area, because it would require much more space. Instead, here we only mention in 
brief the main features of a new recent approach to constructing integrable maps which was 



introduced in g, |T0|, |T2|, developed in 0, which has also been used in this paper and which 
will be referred to as Backlund transformations for finite-dimensional integrable systems. 

One of the new features of this approach to disrete-time integrability is the spectrality 
property which is a projection on the classical case of the famous (quantum) Baxter equation. 
It was discovered on the examples of Toda lattice and elliptic Ruijsenaars system in ^ and 
was generalized to the integrable case of the DST model in |jlO[. Later on, it was observed 
that the property is universal and that, in effect, it gives a canonical way of parameterising 
the corresponding shift on the Jacobian which is characterized by adding a point (A, fi) to a 
divisor of points on the spectral curve F (cf. 0). 

A direct consequence of the spectrality property is the explicitness of the constructed 
maps. This new point, which is an obvious advantage because explicit iterative maps are 
much more useful than implicit maps (given as a system of non-linear equations), was clearly 
demonstrated in 0. This new aspect of constructing explicitly given maps has been also 
adopted and illustrated in detail in the present paper. 

There were several examples of explicit maps known before, like McMillan's map, but all 
those cases were exceptional, for in the generic situation, according to Veselov's approach, 
integrable Lagrange/Poisson correspondences are multi- valued maps, i.e. correspondences 
rather than maps. Using the spectrality property as extra data allows to overcome this 
drawback and to construct discrete-time integrable flows as genuine maps. 

Another new feature of the proposed construction of integrable time-discretizations is an 
identification of the most elementary, one-point, basic map and construction of composite 
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maps, like the two-point map, as compositions of the one-point map and its inverse. The 
choice of the matrix M{u) ( p.8|) generating the one-point map is dictated by the algebraic 
considerations explained in [jl2|. In brief, the matrix M{u) should be a simple L-operator of 
the quadratic algebra, 

{L^{u),L2iv)} = [r{u-v),L,{u)L2iv)], (8.1) 

with the same rational r-matrix ( |2.7D as in the linear algebra ( [^.6[ ). The number of zeros of 
the det M{u) is the number of essential Backlund parameters, so that the matrix M{u) in 
( |3.8D is one-point and the matrix M{u) in ( [7.22[ ) is two-point. The fact that the right ansatz 
for the matrix M{u) obeys the algebra ( pTTI) usually garanties that the resulted map will be 
Poisson, see [|12| for details. 



In the present paper we have observed a new 'spectrality' property of the basic one-point 
map with respect to the parameter 7 in 

detM(M) = 7(m- A). (8.2) 

We have also shown that the two-point map factorises to two one-point maps. 

The two-point map constructed above is probably most general map for the considered 
sl(2) Gaudin model, meaning that it gives a discretization of continuous flows given by any 
Hamiltonian Hj, j = 1, . . . ,n, from the spectral curve, 

= A\u) + B{u)C{u) =a' + J2 f-^ + ^^-^] . (8.3) 

So, at least in principle, any other integrable map for this model should be a function of the 
n maps constructed in this paper. 

There is no established name for integrable maps with all the qualities mentioned above, 
namely: i) spectrality, ii) explicitness, iii) Poissonicity, iv) limits to continuous flows, v) 
preservation of the same integrals as for the continuous flows which these maps discretize. We 
are using for them the same name, Backlund transformations, as was used in the references 

[i, Sim- 

The application of the constructed maps as exact numerical integrators of the continuous 
flows is considered in iK 
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